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Two-loop helicity amplitudes for fermion-fermion scattering 
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We discuss the calculation of two-loop helicity amplitudes for quark-quark scattering in QCD and four-gluino 
scattering in M = 1 supersymmetric Yang-Mills theory. We study the dependence of the results on different 
variants of dimensional regularization. In particular, we consider the 't Hooft-Veltman and four-dimensional 
helicity (FDH) schemes. We also discuss ambiguities in continuing the Dirac algebra to D dimensions. For the 
four-gluino case, once the infrared divergent part of the amplitude is subtracted, the finite remainder is free of 
these ambiguities. 



The level of precision that will be achieved in 
experiments at the LHC poses a challenge for the- 
orists to match. In principle, new physics (or 
hints of it) may be found in small discrepancies 
between theoretical predictions and experimen- 
tal results and it is therefore important that the 
theoretical predictions be robust. As one com- 
ponent of this, it would be helpful to know the 
cross section for the production of hadronic jets 
through next-to-next-to- leading order (NNLO) in 
the QCD coupling. There are a number of im- 
provements to be expected from such calcula- 
tions, for example, reducing the renormalization 
and factorization scale uncertainties in produc- 
tion rates. It should also allow a better under- 
standing of energy flows within jets, as a jet may 
consist of up to to three partons at NNLO. It 
also allows for better matching between parton- 
level and hadron-level jet algorithms. Very im- 
portantly, it is also the first order where hon- 
est assessments of theoretical uncertainties can 
be made. See, for example, ref. pQ for further 
discussions. 

There has been great progress in the past few 
years in the calculation of two-loop matrix el- 
ements, especially for 2 — > 2 scattering pro- 
cesses 0-^)]. This progress has been possible 
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thanks to new developments in loop integration 
|17)-|23j and in understanding the infrared diver- 
gences of the theory [231-1221 ■ br particular, in 
a series of papers the Durham group calculated 
the interference of the two-loop amplitudes with 
the tree level ones, summed over all external color 
and helicity states, for all 2 — * 2 parton processes 
0-0 

Two-loop amplitudes have also been calculated 
keeping full information on the color and helicity 
states of the external particles CDJ- CHI ■ This ad- 
ditional information is not needed for the main 
phcnomenological application, namely, NNLO jet 
production in collisions of unpolarized hadrons. 
However, experiments at RHIC involve the scat- 
tering of polarized protons, for which the helicity 
amplitudes are directly relevant. Other advan- 
tages of having the amplitudes in a helicity basis 
include the study of a number of formal properties 
of scattering amplitudes, such as supersymmetry 
Ward identities |27]-[2E|) collinear and high en- 
ergy behavior (2H-G2], an d the link between color 
decomposed QCD amplitudes, twistor space and 
topological string theories, recently uncovered by 
Witten |33] . 

We present here a summary of our calculation 
for quark-quark scattering amplitudes in QCD, 
as well as four-gluino scattering in M — 1 super- 
Yang-Mills theory |33|. Our results for the two- 
loop quark-quark helicity amplitudes agree with 
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Glover's results ^H] (after correction of minor er- 
rors). 

The three QCD processes considered here are, 
?(Pi,Ai) + g(p 2 ,A 2 ) ^0(P3,A 3 ) + Q(p 4 ,A 4 ), (1) 
g(p l ,Ai) + g(p2,A2)->g(pa,A 3 ) + Q(p4,A4), (2) 
g(p 1 ,Ai) + Q(pa,A2)->g(p3,A 3 ) + Q(p4,A4). (3) 

For the four-gluino case we consider 

g(pi, Ai) + g(p2, A 2 ) -> ff(ps, A 3 ) + jr(p 4 , A 4 ), (4) 

where we use a "standard" (not all outgoing) con- 
vention for the external momentum (p,) and he- 
licity labeling (Ai). 

The amplitudes are calculated using dimen- 
sional regularization. We use the following pre- 
scription when a trace of the Minkowski metric is 
encountered, 

V % = D a =4-2e6 R . (5) 

In this way, we have a continuous set of schemes 
labeled by Sr. Setting Sr = 1 corresponds to 
the 't Hooft-Veltman scheme [35], while setting 
5r = corresponds to the FDH scheme |36I28| . 
The FDH scheme has improved supersymmetry 
properties by virtue of fixing the number of inter- 
nal gluon states to two. 

We evaluate the amplitudes using the spinor 
helicity formalism |37) . In this formalism, 
the amplitudes will be proportional to helicity- 
dependent phase-containing factors written in 
terms of spinor inner products. These spinor in- 
ner products are defined as (ij) = and 
[y] = where \i^) are massless Weyl 

spinors of momentum pi, labeled with the sign 
of the helicity. 

The L-loop amplitudes are color decomposed 
as, 

3 

M {L) =SxY, Tr [c] x jtfW'M , (6) 

c=l 

where S is the spinor factor mentioned above. For 
example, for process Q, with Ai = A4 = + and 
A 2 = A 3 = -, we have S = i(31)/(42). The 
factors for other helicity configurations and for 
process Q are similar. The quantities M^'W 
depend only on the Mandelstam variables s = 
(Pi +P2) 2 , t = (pi - Pi) 2 , and u = {pi -p 3 ) 2 . 



For the quark process Q the color basis is 

TrW =S\6%, T,M=5%5%. (7) 

The color bases for processes (0 and are sim- 
ilar. For the four-gluino case the color decom- 
position is the same as for the four-gluon case, 
and is given in terms of traces of color matri- 
ces (or products of them). Using the notation 

tT (Ta iT a jT a kT aij = we can write 

Tr [1] = tr 1234 , TrW = t r 1243 , Tr^ = tr 1423 

Tr [4] = tr 1324 , Tr [5] = tr 1342 , Tr [6] = tr 1432 , 

TrW = tr 12 tr 34 , Tr® = tri3 tr 24 , 

Tr® = tr 14 tr 23 . (8) 

The two-loop Feynman diagrams were gener- 
ated using QGRAF 38 . A MAPLE program was then 
used to evaluate each diagram. Some of the di- 
agrams were also evaluated using FORM (21) as a 
cross-check. 

When the interference method is used, i.e., 
when one calculates the two-loop amplitudes in- 
terfered with the tree-level ones, summed over 
helicities and colors, one can use standard trace 
techniques to put the integrand into a form con- 
taining only dot products of momenta. Integral 
reduction algorithms then give the loop integrals 
in terms of a minimal set of master integrals. In 
our case we want to keep full information over 
helicity and color states. In order to put the in- 
tegrals into a form suitable for applying the gen- 
eral reduction algorithms, we multiply and divide 
by appropriate spinor inner products constructed 
from the external momenta. These spinor in- 
ner products effectively play the role of the tree- 
level amplitudes in the interference method, ex- 
cept that in this method the helicity information 
is maintained when converting the spinor strings 
into traces over 7 matrices. 

An important way of checking the correctness 
of the calculation is by comparing the infrared 
divergence of the renormalized amplitudes with 
the ones predicted by Catani's formula for two- 
loop n-point amplitudes 24 , 

|^(2)) =J W | M W) + J(2) |^(0)) 

+ lM 2)fin ), (9) 
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where the "ket" notation \A4^} indicates that 
the L-loop amplitude is treated as a vector in 
color space. The components of this vector are 
given by the M^'^°' of eq. ©. The divergences 
of M.^ are encoded in the color operator 

(2) 

while those of Aih, also involve the scheme- 
dependent operator j' 2 '. Catani's formula 10 
not only allows us to check the results, but also 
to organize them in terms of finite and divergent 
parts. Proofs of Catani's formulas have appeared 
in refs. \£A%\ . 

For each process and each color basis we will 
have a different 1^ matrix. For the basis we 
have, 



= e-<+W ( 1 :! 



r(l-e) \e 2 2e 

T - U 

2C F S-^(T-U) J ' 



(10) 



2cyr-^(s-u 
s -u 



where N is the number of colors, Cf = (N 2 
1)/(2N) and 



U = 



The corresponding operator for qQ qQ is ob- 
tained by changing S — > U, T — » S and U — ► T in 
eq. Similarly, the operator for qQ — » qQ is 

obtained by exchanging S and U in (|10|) . For the 
gluino case we get a nine-by-nine matrix, given in 
eq. (2.18) of ref. [Ej. 

The operator 1^ is given by |23] 



r(2) 



e+^)r(l-2e) (bp } (1) 



(2) 



(11) 



where bo is the first coefficient of the QCD /3- 
function, and the coefficient K depends on 5r 
and is given by |24I13| 
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+ 7^(1-**) 



(12) 



with Afy- being the number of massless fundamen- 
tal representation quarks. is a universal op- 
erator of order 1/e. A general expression for J? < - 2 - ) 
with an arbitrary number of external legs was re- 
cently presented in |2(i| . 

The two-loop remainders, as defined by eq. 10 
will have the following form 



2), [elfin 



Qx-K\ In 
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M (o), 



+ Q 2 



(l)i[c]fin 



,[c] 



(13) 



where Qi, Q2 and Q3 depend on the number 
of colors N, number of flavors Nf, and scheme 
label 5r (in particular, Q 1 = Q 2 = Q3 = 
when 5r = 1). The parameter fj, is the renor- 
malization scale. The M^ 1 )'^ 11 are the one-loop 
finite remainders, and the M( 1 ' , >- C ' ,€ ' s * t are the 
^-dependent parts of the 0(e 1 ) coefficients of 
the one- loop amplitudes. Finally, consists of 
powers of N and Nf multiplying functions (loga- 
rithms, poly logarithms, as well as rational func- 
tions) of the Mandelstam variables. The explicit 
form of these functions is given in ref. |34|. 

For the four-gluino case in Af = 1 super- Yang- 
Mills theory, the two-loop finite remainders in the 
FDH scheme (Sr = 0) are, 



M (2 



-6g 



In 



" K In 



r(0),[c] 



for c = 1 



-2b (ln^)-in)M^ 

+N 2 AW + , (14) 
. 6, and 

/7 rh/ l "' [c]fin 

(15) 



-2b \ \a 
f NG [c] , 
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for c = 7,8,9. Here bo and b\ are the first two 
coefficients of the J\f = 1 super- Yang-Mills (3- 
function. AM, and are functions of s, 
t and it. Again the explicit form of all these func- 
tion is given in ref. |34| . 

In previous papers |28I13I14] we showed that 
the following supersymmetry Ward identities are 
satisfied at two loops through 0(e°) when one 
uses the FDH scheme: 



(16) 
(17) 



M slJSY (gf. 


■92 


,9$ ■ 


•9i 


M slJSY (gt. 


92 


,9a: 


si 


M SUSY (gt. 


92 


>9s ■ 


gt 



(23) 
(13) 



M 



SUSY 



{gt,9-2 ,9s >9t)' (18) 



With the four-gluino amplitudes, we may also 
check the supersymmetry identity relating the 
four-gluon amplitude to the four-gluino one, 

^ SY (9t, 92,9^-91) 
(2 4) 



(13) 



A4 SUSY ( 5 + s^,^). (19) 



(Note that in ref. ^3] an aU outgoing definition 
of helicity is used for the four-gluon amplitudes, 
while we use the 'standard' one here where legs 
1 and 2 are incoming.) It turns out that this 
last identity does not work immediately at two 
loops |34|. The problem is related to ambiguities 
arising from continuing the Dirac algebra to D 
dimensions. At one-loop the ambiguity is harm- 
less because it affects only 0(e) terms. From 
the Catani formula a one-loop ambiguity at 0(e) 
leads to an ambiguity in the 0(l/e) terms at two 
loops. It is of course possible to arrange for a 
prescription to fix the ambiguities to restore the 
manifest supersymmetry Ward identities, but, in 
any case, these two-loop ambiguities may all be 
absorbed into Catani's formula for the divergent 
parts, leaving well defined finite parts. As ex- 
pected, in the FDH scheme the finite remainders 
A, B and G in eqs. (|f4|l and l|15|l satisfy supersym- 
metry identities and agree with the corresponding 
functions given in ref. |13j for pure glue scatter- 
ing. 
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